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Field of the Invention 

This invention relates generally to image processing and, more specifically, to 
15 detection of cloud cover in high- altitude and/or orbital overhead imaging data. 

Background of the Invention 

Overhead imaging studies of a surface below may be hampered by the presence of 
cloud formations. Understandably, thick clouds between an observation point and the area of 
interest under observation can conceal objects or features in the area of interest. Potentially 

20 worse in some cases is the presence of thinner cloud formations that do not entirely occlude 
the surface, but may reduce the contrast of surface features and change the derived surface 
spectral reflectance signature with resulting impact on information products such as spectral 
vegetation indices. Presence of thin cloud formations, such as cirrus clouds, can skew the 
analysis of such surface features by causing researchers to confuse presence of cloud features 

25 for features or changes in the surface region of study. For example, FIGURE 1A shows a 
representative image 100 of a surface area under study. Merely looking at the image, it may 
be difficult to determine which aspects of the image are surface features 110 and which 
aspects are cloud features 120. 

Because the presence of cloud formations can interfere with the accuracy of overhead 

30 imaging studies, methodologies have been developed to detect the presence of cloud 
formations so that accuracy of surface studies will not be undermined by undetected cloud 
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patterns. One approach is to use "clear-sky" spectral or reflectance maps of the areas of 
interest to detect the presence of clouds. By comparing the clear-sky maps with current 
imaging data, large-area spectral or reflectance changes may signal the presence of cloud 
cover. This approach involves successfully collecting, verified clear-sky imaging data of the 
5 area of interest. The clear-sky maps typically are created using thermal infra-red 
measurements to determine the presence of cloud formations. Most cloud formations, 
including high altitude cirrus clouds made up of ice crystals, present a distinct, differentiable 
thermal signature. If thermal data indicates the presence of cirrus or other clouds in an area 
of study, it will be understood which portions of the image data are affected by the presence 

10 of clouds. Thus, analysis of the area of interest will not be distorted by the presence of 
undetected cloud formations. 

FIGURE IB shows a "cloud mask" 150 derived using conventional techniques to 
show the cloud features 120 in the original image 100 of FIGURE 1A. Absent the cloud 
mask 150, it can be appreciated that it might have been easy to confuse edges of cloud 

1 5 patterns 1 20 with surface features 110. 

Unfortunately, as is readily appreciated, collection of thermal-infra red data requires 
equipment capable of gathering thermal-infrared data. In the realm of orbital satellites, 
integrating such equipment into the satellite increases cost. Additional telemetry involved in 
making use of such data also is resource-intensive and costly. 

20 Even where such clear-sky data are available, continual accurate analytical 

comparison of archival clear-sky data with currently-captured imaging data is needed to 
ensure that the captured data represents suitably accurate, cloud-free images. Determination 
of whether the imaging data is suitably cloud-free is a significant concern. If it is not 
accurately determined whether captured images are suitably cloud-free, it may be necessary 

25 to arrange for the areas of interest to be re-imaged. Analysts who desire to use images from 
an image archive need to be assured that the image data is sufficiently cloud-free to be 
worthy of acquisition and use in their research. In addition, before quantitative analysis tools 
are applied to analyze the imaging data, the imaging data must be determined to be suitably 
cloud- free to ensure that the resulting quantitative analyses will be correct. 

30 Thus, there is an unmet need in the art for a method for determining presence of 

clouds in aerial imaging data not involving use of special thermal infra-red sensing 
equipment or the data collected by such equipment. 

Summary of the Invention 

Embodiments of the present invention can be used to determine the presence of 
35 clouds without involving thermal infrared detecting equipment. Embodiments of the present 
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invention use spectral and spatial tests applied to pixel-level spectral measurements to 
determine the presence of cloud formations. The tests are computationally simple and, thus, 
do not impose an unreasonable operational computing workload. A sequence of such tests 
may be successively applied to the pixel-level spectral measurements to classify the pixel as 
5 indicating presence or absence of a cloud. 

The present invention comprises methods, a computer-readable medium storing 
instructions, and a system for determining whether a data point of an image indicates a 
presence of a cloud using data including visible, near-infrared (NIR), and short wavelength 
infrared (SWIR) data. In one embodiment, a first comparison of a cirrus-band reflectance of 

10 a data point with a threshold cirrus-band reflectance value is made, classifying the data point 
as a cloud point if the cirrus-band reflectance of the data point exceeds the threshold cirrus- 
band reflectance value. When the comparing of the cirrus-band reflectance of the data point 
with the threshold cirrus-band reflectance value does not sufficiently classify the data point 
as a cloud point, a further analysis of the data point is performed. The further analysis 

1 5 includes performing a second comparison of an additional cloud indicator with an additional 
cloud indicator threshold, the additional cloud indicator being derived from at least one of the 
visible, near-infrared, and/or short wavelength infrared data. The data point is classified as 
one of a cloud point or a non-cloud-point when the second comparison of the additional 
cloud indicator with the additional cloud indicator threshold allows the data point to be 

20 classified as one of a cloud point or a non-cloud point. 

Brief Description of the Drawings 

The preferred and alternative embodiments of the present invention are described in 
detail below with reference to the following drawings. 

FIGURE 1A is a conventional aerial image of an exemplary region of interest in 
25 accordance with the prior art; 

FIGURE IB is a conventional cloud mask derived from conventional techniques to 
indicate the presence of clouds in the image of the exemplary region of interest in accordance 
with the prior art; 

FIGURE 2 is a graph plotting NDSI values versus D values and illustrating points 
30 where clouds are present; 

FIGURE 3 is a flowchart of a routine according to an embodiment of the present 
invention for determining presence of cloud formations; 

FIGURE 4 is a flowchart of a routine according to another embodiment of the present 
invention for determining presence of cloud formations; and 
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FIGURE 5 is a block diagram of a system according to an embodiment of the present 
invention. 

Detailed Description of the Invention 

By way of overview, the present invention comprises methods, a computer-readable 
5 medium storing instructions, and a system for determining whether a data point of an image 
indicates a presence of a cloud using data including visible, near-infrared, and short 
wavelength infrared data. In one embodiment, a first comparison of a cirrus-band reflectance 
of a data point with a threshold cirrus-band reflectance value is made, classifying the data 
point as a cloud point if the cirrus-band reflectance of the data point exceeds the threshold 

10 cirrus-band reflectance value. When the comparing of the cirrus-band reflectance of the data 
point with the threshold cirrus-band reflectance value does not sufficiently classify the data 
point as a cloud point, a further analysis of the data point is performed. The further analysis 
includes performing a second comparison of an additional cloud indicator with an additional 
cloud indicator threshold, the additional cloud indicator being derived from at least one of the 

15 visible, near-infrared, and/or short wavelength infrared data. The data point is classified as 
one of a cloud point or a non-cloud-point when the second comparison of the additional 
cloud indicator with the additional cloud indicator threshold allows the data point to be 
classified as one of a cloud point or a non-cloud point. 

Studying images, empirically it can be determined for each of these data points 

20 whether the data point signifies a cloud point or a non-cloud point. It will be appreciated 
that, in accordance with embodiments of the present invention, a number of quantities can be 
calculated for each data point using data extractable from visible, near-infrared, and short- 
wavelength infrared data. By studying these calculated quantities, threshold values are 
determinable by which the calculated quantities suitably are used to automatically determine 

25 whether a data point represents a cloud point or a non-cloud point. It will also be appreciated 
that, although embodiments of the present invention may analyze data to determine a 
presence of both visible clouds and sub-visible cloud layers, a presently preferred 
embodiment of the present invention is tailored to classifying data points based on whether 
the data points indicate the presence of visible clouds. 

30 FIGURE 2 shows an exemplary graph 200 of such calculated quantities. Examining 

such a graph or other data representation in conjunction with associated cloud masks and 
RGB image, it can be determined what data points should be classified as cloud points and 
which should be classified as non-cloud points. More specifically, the graph shows 
normalized difference snow index, NDSI, plotted on a vertical axis 210 against values of a D 
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variable on the horizontal axis 220. In one embodiment, NDSI is determined by equation 
(1): 

NDSI = (p Green -PswiRx)l{pGreen+PsWlRx) 0) 

5 

The reflectance values, PGreen and Pswiri, represent reflectance in the selected wavelength 
range, such as the green wavelengths, the short-wave infrared wavelength, respectively. The 
D variable, in turn, is determined from a normalized difference vegetation index, NDVI, 
respectively determined from equations (2) and (3): 

10 

D = \NDVl\°- 6 /(p K J (2) 
ND VI = (p N/R - p Red )/{p NIR + p Re , ) (3) 

15 The reflectance values, pR e d and p N i R , represent reflectance in the selected wavelength range, 
such as the red wavelengths, the near-infrared wavelengths, respectively. 

The graph 200 shows values of NDSI versus D for data points that have been 
empirically identified by manual, visual inspection of an area of interest as representing 
bodies of water 230, ground 250, and clouds 260. Thus, because the data used in deriving 

20 NDSI, D, and NDVI is derivable from visible/near-infrared and short-wavelength infrared 
data, data points representing cloud points and non-cloud points can be identified without 
separate thermal infrared data. 

A comparable analysis is achievable using different formulations of D. For example, 
in analyzing data collected by the Multiangle Imaging SpectroRadiometer (MISR) sensor 

25 used by NASA, D is calculated by raising NDVI to different exponential values depending 
on a type of ground cover expected to be present in the surface being imaged. Although the 
MISR D use is more complex because of its landcover-type-dependent NDVI exponent and 
large, statistically derived, D-threshold database, MISR D values also can be used with 
embodiments of the present invention to achieve satisfactory results. 

30 Embodiments of the present invention can employ a number of such quantities to 

classify data points as cloud points or non-cloud points. Selection, ordering, calculation, and 
comparison of such quantities can be made in order to balance computational burdens and 
desired classification precision. For example, in two exemplary embodiments described 
below, a first comparison involves reflectance in the cirrus-band wavelengths, p C u with a 

35 threshold value which provides a ready, reliable first step in classifying data points as either 
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cloud points or non-cloud points. In contrast with NDSI or D, pci can be compared to a 
threshold value without additional computation, thus making a comparison of pa as a first 
step can reduce computational burdens. It will also be appreciated that the steps can be 
ordered to evaluate potentially more reliable classifiers first, or the steps can be ordered to 
5 provide a logical AND or OR construct to provide for reliable classification of the data 
points. 

Depending upon the computational resources available, it will be appreciated that 
analysis of data points can occur in real-time, while analysis of classification precision vs. 
computing load may occur in non-real-time.. If non-real-time analyses indicate that greater 

1 0 precision can be achieved, to better meet user needs within available computing resources, by 
adjusting thresholds or adding additional cited tests to the real-time test hierarchy for specific 
background landcover types, locations, or times of year, those revisions can be made for 
future real-time processing. 

FIGURE 3 shows a first embodiment of a routine 300 for classifying data points as 

15 cloud points or non-cloud points. The routine 300 begins at a block 302, and a next point of 
top of atmosphere, visible, near-infrared, and short-wavelength infrared data is submitted for 
processing at a block 304. At a block 306, a comparison of the p C i to a cirrus-band threshold 
cloud value is made. If pci exceeds the threshold value, the data point is classified as a cloud 
point at a block 308. If not, the routine 300 proceeds to a next block to make a further 

20 comparison of whether the data point represents a cloud point or a non-cloud point. 

In one particular embodiment, the p C i comparison at the block 306 is made at a 
wavelength of 1.88 jim. At this wavelength, the reflectance has been determined to be more 
reliable than at slightly lower wavelengths. Of course, in alternate embodiments, p C i may be 
tested at wavelengths other than 1.88 //m, such as at 1.38 fim or other cirrus bands. 

25 If the comparison of p C i at the block 306 to make an initial determination of whether 

the data point was a cloud point did not result in the data point being classified as a cloud 
point at the block 308, additional comparisons can be made to further differentiate whether 
the data point is a cloud point or a non-cloud point. The comparisons and number of 
comparisons selected suitably are chosen to balance between computational simplicity and 

30 classification precision. Generally, as a greater number of comparisons are performed, 
greater precision is obtained. Nonetheless, selecting a fewer number of comparisons may 
result in a desirable degree of accuracy with fewer comparisons and/or calculations being 
made. 

More specifically, if the comparison of the data point at the block 306 does not result 
35 in the data point being identified as a cloud point at the block 308, at a block 3 10 the NDSI is 
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compared to an NDSI snow threshold value. This comparison may eliminate data points 
showing snow. If the data point NDSI is greater than the NDSI snow threshold value, the 
data point is a snow point. Again, the NDSI threshold value may be empirically determined 
using other information from which data points have previously been classified as cloud 
5 points or non-cloud points. If at the block 310 the NDSI exceeds the NDSI snow threshold 
value, the data point is classified as a non-cloud point at the block 312. 

It will be appreciated how threshold values like the NDSI snow threshold value 
compared at the block 310 can affect classification precision. If, for example, the NDSI 
snow threshold is lowered, more data points may be classified as non-cloud ground points. If 

10 analysis reveals that this revision results in a net improvement in classification accuracy, 
application of further comparisons in the routine 300 may be avoided. Adjusting the 
thresholds in the tests described will determine how the individual data points in thinly cloud 
covered areas are classified. Accordingly, selection of thresholds based on empirical analysis 
of tested values for test data points known to be cloud points or non-cloud points over 

15 specific categories of landcover, location and season will incorporate a predetermined 
classification accuracy into embodiments of the present invention. 

If the comparison of the data point at the block 310 does not result in the data point 
being identified as a non-cloud point at the block 312, then at a block 314, a comparison of a 
ratio of the near infrared data to the short- wavelength infrared data, NIR/SWIR1, to a 

20 NER/SWIR1 snow threshold value is made to potentially eliminate data points showing snow. 
If the NIR/SWIR1 value exceeds the NER/SWIR1 snow threshold value, the data point is 
classified as a non-cloud point at the block 312. 

If the comparison of the data point at the block 314 does not result in the data point 
being identified as a non-cloud point at the block 312, then at a block 316 a comparison of a 

25 ratio of the NDSI value to an NDSI cloud threshold value is made to potentially eliminate 
data points showing bright ground. If the NDSI value is less than the NDSI cloud threshold 
value, the data point is classified as a non-cloud point at the block 312. 

If the comparison of the data point at the block 316 does not result in the data point 
being identified as a non-cloud point at the block 312, at a block 318 a comparison of the D 

30 variable is made with a D variable cloud threshold to potentially eliminate data points 
showing vegetation . If the D variable is greater than the D variable cloud threshold, the data 
point is classified as a non-cloud point at the block 312. 

If the comparison of the data point at the block 318 does not result in the data point 
being identified as a non-cloud point at the block 312, then at a block 320 a comparison of a 

35 D spatial variability index, DSVI, is made with a DSVI cloud threshold to potentially 
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eliminate data points showing non-smooth features, 
spatial variability index may be given by: 



In one particular embodiment, the D 



DSVI = \D m -D c \ 



(4) 



5 

D m is mean of D values for at least a three-by-three matrix of data points surrounding the 
data point and D c is a central pixel in the three-by-three matrix of data points. If the DSVI is 
greater than the DSVI cloud threshold value, the data point is classified as a non-cloud point 
at the block 312. 

10 It will be appreciated that calculation of the DSVI is a more computationally 

intensive step than other steps previously undertaken. The DSVI is derived from a plurality 
of D values which, in turn, are calculated from reflectance data of the data point. It will be 
appreciated that this step is not a first step in the routine 300 allowing for the possibility of 
faster, less-intensive methods associated with the foregoing blocks allowing for the data 

15 point to be classified as a cloud point 308 or a non-cloud point at the block 312. On the other 
hand, should additional computing power be available, the D m portion of DSVI could be 
computed for a larger matrix of points such as a mean of a five-by-five or larger matrix, 
centered on D c . Use of a larger matrix can increase the accuracy of the DSVI comparison by 
providing a statistically better D m portion. If the comparison of the data point at the block 

20 320 does not result in the data point being identified as a non-cloud point at the block 312, 
then at a block 322 a comparison of the short-wavelength reflectance, Pswiri, is made to a 
short-wavelength reflectance cloud threshold to potentially eliminate data points showing 
dark features. If Pswiri is less than the short-wavelength reflectance cloud threshold, the data 
point is classified as a non-cloud point at the block 312. 

25 If the comparison of the data point at the block 322 does not result in the data point 

being identified as a non-cloud point at the block 312, then at a block 324 a comparison of 
pRed to a red wavelength cloud threshold value is made to eliminate additional data points 
showing dark features. If p Re d is less than the red wavelength cloud threshold value, the data 
point is classified as a non-cloud point at the block 312. 

30 If the comparison of the data point at the block 324 does not result in the data point 

being identified as a non-cloud point at the block 312, then at a block 326, a comparison of a 
ratio of the NIR/SWIR1 to a NIR/SWIR1 cloud threshold value is made to potentially 
eliminate additional data points showing bright ground. If the NIR/SWIR1 value is less than 
the NIR/SWIR1 cloud threshold value the data point is classified as a non-cloud point at the 

35 block 312. On the other hand, if the NIR/S WIR1 value is greater than the NIR/SWIR1 cloud 
threshold value, the data point is classified as a cloud point at the block 308. 
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Once the data points have been classified as one of a cloud point at the block 308, or 
as a non-cloud point at the block 312, then at a block 328 it is determined if all data points of 
interest have been classified. If not, the routine 300 loops to the block 304 where the next 
data point is addressed. However, if it is determined at the block 328 that all the data points 
5 of interest have been analyzed, the routine 300 ends at the block 330. 

As previously described, the routine 300 uses cloud thresholds empirically derived 
from manual or other studies of overhead imaging data. The threshold values may vary 
depending on the nature of the area of interest and the season during which the imaging data 
is captured. For example, threshold values for forests or closed shrub areas will vary 
10 between summer/tropical seasons and snowy seasons, just as the threshold values will vary 
between permanent wetlands and permanently snow-covered areas. For example, and not by 
way of limitation, Table 1 presented below lists representative threshold values that suitably 
may be used in the routine 300 for scenes and seasons of interest: 



Scene 


Season 


D 


NDSI snow 


NDSI c!ood 


DSVI 


PSWIRl 


NIR/SWIRl snow 


NIR/SWIRl cloud 


PRed 


Pcioud 


Forest/ 
Closed 
Shrub 


Summer 
or 

Tropical 


20 


N/A 


-0.25 


1 


0.1 


N/A 




0.1 


0.03 




Spring or 
Fall 

without 
snow 


10 


N/A 


-0.35 


3.5 


0.1 


N/A 




0.1 


0.03 


««» 


Fall, 

Spring or 

Winter 

(snow) 


2 


0.55 


-0.3 


0.2 


0..2 


4 


N/A 


0.1 


0.03 


Grass or 

Crops 

(Mosaic) 


Summer 


20 
(40) 


N/A 


-0.3 


1 


0.1 


N/A 


1 


0.1 


0.03 




Spring or 
Fall w/o 
snow 


10 


N/A 


-0.35 


1 


0,1 


N/A 


0.9 


0.1 


0.03 




Fall, 

Spring or 

Winder 

w/snow 


2 


0.55 


-0.35 


0.2 


0.2 


4 


1 


0.2 


0.03 


Snow 
and Ice 


All 


2 


0.55 


-0.3 


0.2 


0.2 


4 


N/A 


0.1 


0.03 


Barren 
or Sparse 
Open 
Shrub 


(No 
Snow) 


4 


N/A 


-0.25 


1 


0.1 


N/A 


0.8 


0.2 


0.03 


Savanna 


(No 
Snow) 


10 


N/A 


-0.35 


1 


0.1 


N/A 


0.8 


0.2 


0.03 



15 Table 1 



It will be appreciated that thresholds can be derived from study of other scenes and 
terrains, such as wetlands or water-covered areas as well. 

Various combinations of tests can be used to optimally balance desires for accuracy 
20 and computing efficiency. For example, accurate results are obtainable according to a subset 
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of the routine 300 (FIGURE 3) where comparisons at the decision blocks 306, 316, 318, 320, 
and 326 and at least one of the comparisons at the decision blocks 310, 314, and 322. Table 
2, on a next page, shows a computed accuracy for tests and combinations of tests run on a 
number of data sets. 
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FIGURE 4 shows a second embodiment of a routine 400 for classifying data points as 
cloud points or non-cloud points. The routine 400 begins at a block 402, and a next point of 
top of atmosphere visible, near-infrared, and short-wavelength infrared data is submitted for 
5 processing at a block 404. At a block 406, a comparison of the pci to a cirrus-band threshold 
cloud value is made. If pci exceeds the cloud threshold value, the data point is classified as a 
cloud point at a block 408. If not, the routine 400 proceeds to a next block to make a further 
comparison of whether the data point represents a cloud point or a non-cloud point. 

If the comparison of the data point at the block 406 does not result in the data point 
10 being identified as a cloud point at the block 408, then at a block 410 the NDSI is compared 
to an NDSI minimum threshold value and an NDSI maximum threshold value. The NDSI 
thresholds are set according to empirical analysis of NDSI data such as that shown in 
FIGURE 2. Unlike the routine shown in FIGURE 3 which uses single-value cutoff 
thresholds, at the block 410 the threshold defines an area of the graph 200 (FIGURE 2) as 
15 opposed to an intercept. More specifically, the comparison of NDSI at the block 410 is: 

NDSI > [(Nm-N l )/D t ]*D+N l 
AND 

NDSI < [(Nm-Nh)/D t ]*D+N l 

20 

If either comparison is false, the data point is classified as a non-cloud point at a block 412. 

Values for these threshold calculation numbers are included in Table 2, below. 

If the comparison of the data point at the block 410 does not result in the data point 

being identified as a non-cloud point at the block 412, at a block 414 a comparison of a ratio 
25 of NIR/SWIR1 to a NIR/SWIR1 snow threshold value is made. If the NIR/SWIR1 value is 

greater than the NIR/SWIR1 snow threshold value, the data point is classified as a non-cloud 

point at the block 412. 

If the comparison of the data point at the block 414 does not result in the data point 

being identified as a non-cloud point at the block 412, then at a block 416 a comparison of a 
30 ratio of the NIR/SWIR1 value to an NIR/SWIR1 cloud threshold value is made. If the 

NIR/SWIR1 value is less than the NIR/SWIR1 cloud threshold value, the data point is 

classified as a non-cloud point at the block 412. 

If the comparison of the data point at the block 416 does not result in the data point 

being identified as a non-cloud point at the block 412, then at a block 418 a comparison of 
35 p R ed to a red wavelength cloud threshold value is made. If p Re d is less than the red 
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wavelength cloud threshold value, the data point is classified as a non-cloud point at the 
block 412. 

If the comparison of the data point at the block 418 does not result in the data point 
being identified as a non-cloud point at the block 412, then at a block 420 a comparison of 
5 the short-wavelength reflectance, Pswiri, is made to a short-wavelength reflectance cloud 
threshold. If Pswiri is less than the short-wavelength reflectance cloud threshold, the data 
point is classified as a non-cloud point at the block 412. 

If the comparison of the data point at the block 420 does not result in the data point 
being identified as a non-cloud point at the block 412, then at a block 422 a comparison of 
10 the DSVI is made with a DSVI cloud threshold. If the DSVI exceeds the DSVI cloud 
threshold, the data point is classified as a non-cloud point at the block 412. On the other 
hand, if the DSVI is less than the DSVI cloud threshold, the data point is classified as a cloud 
point at the block 408. 

Once the data points have been classified as one of a cloud point at the block 408 or 
15 as a non-cloud point at the block 412, then at a block 424 it is determined if all data points of 
interest have been classified. If not, the routine 400 loops to the block 404 where the next 
data point is addressed. However, if it is determined at the block 424 that all the data points 
of interest have been analyzed, the routine 400 ends at the block 426. 

As previously described, the routine 400 uses threshold calculations empirically 
20 derived from manual or other studies of overhead imaging data. The threshold values may 
vary depending on the nature of the area of interest and the season during which the imaging 
data is captured. For example, threshold values for forests or closed shrub areas will vary 
between summer/tropical seasons and snowy seasons, just as the threshold values will vary 
between permanent wetlands and permanently snow-covered areas. Again, for example and 
25 not by way of limitation, Table 3 below lists representative threshold values that suitably 
may be used in the routine 400: 



Scene 


T(CI) c |oud 


N L 


N M 


N„ 


D T 


T(NIR/SWIRl) snow 


T(NIR/SWIRl) clood 


pRed 


SWIR1 


DSVI 


Forest - 
summer 
mid- 
latitude 


0.03 


-0.2 


-0.15 


1 


20 


N/A 


0.8 


0.1 


0.1 


1 


Forest - 
tropical 


0.03 


-0.5 


-0.3 


1 


20 


N/A 


0.8 


0.1 


0.1 


1 


Crops or 
Mosaic - 
Summer 


0.03 


-0.3 


-0.2 


1 


40 


N/A 


0.8 


0.1 


0.1 


1 


Grass - 
Summer 


0.03 


-0.3 


-0.2 


1 


20 


N/A 


0.8 


0.1 


0.1 


1 


Crops or 
Mosaic - 
Spring 


0.03 


-0.3 


-0.2 


1 


10 


N/A 


0.8 


0.1 


0.1 


1 
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Crops or 
Mosaic - 
Snow 


0.03 


-0.5 


-0.3 


0.6 


5 


4 


0.8 


0.1 


0.1 


1 


Barren 


0.03 


-0.3 


-0.2 


1 


3 


N/A 


0.8 


0.1 


0.1 


1 


Savanna 
or Open 
Shrub 


0.03 


-0.3 


-0.2 


1 


10 


N/A 


0.8 


0.1 


0.1 


1 



Table 3 



It will be appreciated that the routine 400 (FIGURE 4) simplifies the selection of threshold 
5 values. 

FIGURE 5 shows a system 500 according to an embodiment of the present invention. 
Imaging data 5 10, including imaging data from at least one data point, is received. Threshold 
data 520, such as the parameters previously described in connection with FIGURES 3 and 4 
and Tables 1 and 2, is supplied to the system for comparison. A cirrus band comparator 530 

10 makes a first comparison of the data point with a cirrus band threshold. As previously 
described, if the cirrus band reflectance of the data point exceeds the cirrus band threshold, 
the data point is classified as a cloud point in cloud mask data 550 and/or a cloud mask 560. 
On the other hand, if use of the cirrus band comparator 530 does not result in classification of 
the data point, a secondary comparator 540 is applied to classify the data point. Using 

15 routines previously described in connection with FIGURES 3 and 4, the secondary 
comparator 540 uses additional cloud indicators and cloud indicator thresholds to classify the 
data points. When the secondary comparator 540 classifies the data point as either a cloud 
point or a non-cloud point, the data point is appropriately classified in the cloud mask data 
550 and/or the cloud mask 560 in accordance with the predetermined classification precision 

20 determined by the threshold levels established for the comparisons being made. The system 
500 suitably is applied to all data points in the imaging data 510 to generate cloud mask data 
550 and/or a cloud mask 560 for the imaging data 510. 

It will be appreciated that, in one embodiment of the invention, the determination as 
to acceptability of accuracy provided by a given set of tests and thresholds would be 

25 determined by processing a representative set of imagery off-line, and by evaluating resultant 
cloud masks in comparison with "truth" cloud masks produced by expert analysis. 
Adjustments would be iterated and evaluated for optimization, and production test hierarchy 
and thresholds would then be adjusted for subsequent on-line production runs. 

While alternate and preferred embodiments of the invention have been illustrated and 

30 described, as noted above, many changes can be made without departing from the spirit and 
scope of the invention. Accordingly, the scope of the invention is not limited by the 
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disclosure of the preferred embodiment. Instead, the invention should be determined entirely 
by reference to the claims that follow. 
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